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Using methods made possible by recent advances in photoemission technology, we perform an in- 
depth hne-shape analysis of the angle-resolved photoemission spectra of the electron doped (n-type) 
cuprate superconductor Ndi.85Ceo.i5Cu04. Unlike for the p-type materials, we only observe weak 
mass renormalizations near 50-70 meV. This may be indicative of smaller electron-phonon coupling 
or due to the masking effects of other interactions that make the electron-phonon coupling harder to 
detect. This latter scenario may suggest limitations of the spectral function analysis in extracting 
electronic self-energies when some of the interactions are highly momentum dependent. 

PACS numbers; 79.60.Bm, 73.20.Dx, 74.72.-h 



I. INTRODUCTION 

Angle-resolved photoemission spectroscopy (ARPES) 
is one of the most direct probes of the electronic structure 
of solids. In recent years, with the advent of Scienta an- 
alyzers and high-flux beam lines, a new era of photoemis- 
sion line-shape analysis has been ushered in [1,2]. These 
advances offer unprecedentedly high momentum and en- 
ergy resolution, as well as the ability to perform simulta- 
neous parallel angle scanning and thereby generate direct 
two dimensional (2D) E — k images of ARPES spectral 
functions. Whereas, ARPES data have been displayed 
traditionally as energy distribution curves (EDC's) in 
which the photoemission intensity is plotted as a func- 
tion of energy at specific angles, now they can be repre- 
sented in terms of these detailed 2D intensity plots and in 
analogy with EDC's, complementary momentum distri- 
bution curves (MDC's) can be generated. Such advances 
are making new analysis methods possible, where one 
is able to bring to bear the considerable mathematical 
machinery of many-body physics for intuition and inter- 
pretation. 

Within the sudden approximation, ARPES measures 
the single-particle spectral function A{k,bj) [1-3]. As- 
suming that the material under investigation has well- 
defined electronic excitations [4] , the single-particle spec- 
tral function can be represented compactly in terms of a 
complex self-energy E = E' -f iE" as 



TT [w - - E'(A?, w)]2 + [E"(fc, w)]2 



In momentum regions where E(fc,u;) is a weak func- 
tion of k and the bare dispersion can be linearized as 
— iPp ■ {k — kp) (where iPp is the bare band veloc- 
ity) A{k, uj) at constant oj can be put in the form of a 
Lorentzian centered at w/zJ^ -I- kp — ^' {k , lu) / ifj^ with 
width E"(fc,w)/z;0. 

Fits to MDC's gives one, in principle, not only a mea- 
sure of E"(fc,a;), but also a more accurate parameteri- 
zations of the renormalized dispersion in the case of a 
rapidly changing E'(A:,a;). Such an analysis allows one 
to probe the dominant kinds of electron-electron inter- 
actions as well as for signatures of electronic interaction 
with bosonic modes. Such studies are an area of much 
current focus in the ARPES community with, for in- 
stance, the observation of bosonic effects in the ARPES 
spectra of systems with strong electron-phonon coupling 
such as Be and Mo [6-8]. Recently, similar features have 
been found in the low-energy ARPES spectra of a number 
of high-temperature hole doped cuprate superconductors. 
A mass renormalization, or "kink" , in the dispersion has 
been found ubiquitously in the p-type materials (at w 70 
meV) [9,10]; its origin is a matter of much current debate. 
Some groups have pointed to a many-body electronic 
source [11-13], that is perhaps related to the magnetic 
resonance mode discovered via neutron scattering. How- 
ever, this assignment has been disputed on the grounds 
that the resonant mode may have insufficient spectral 
weight to cause the observed effect with reasonable es- 
timates for its coupling strength [14,15]. Others have 
argued that the kink's presence above Tc and its uni- 
versality among material classes demonstrate a phononic 
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origin and indicates the strong role that lattice effects 
have on the low-energy physics [10]. 

Although, in principle, it is straightforward with the 
above analysis to extract the electronic sclf-cncrgy from 
ARPES spectra, it is not as obvious how to assign parts 
of the self-energy to different origins. For instance, is 
the phenomenon that causes the kink the same one that 
causes the pseudogap in the underdoped materials, or are 
multiple effects playing a role? Moreover, there have even 
been reports that distinguish the kinks near the (tt, 0) and 
(7r/2, 7r/2) positions [13]. Similar manifestations of the 
various phononic, antiferromagnetic, charge order, stripe, 
and structural effects (not to mention the possibility of 
hidden orders [16]) that may exist in the cuprates, makes 
it difficult to extract unambiguous information from a 
spectral function analysis. 

Such problems may be particularly acute in the elec- 
tron doped system Nd2-xCea;Cu04 [17]. In previous 
ARPES work on the "optimally-doped" (x = 0.15) sys- 
tem, high-energy pseudogap effects (i.e. the suppression 
of low-energy spectral weight over an extended frequency 
range) similar to those in the underdoped p-type com- 
pounds were observed [18,19]. However, unlike in the 
p-type materials, these regions of maximum pseudogap 
were not on the Fermi surface (FS) near (tt, 0) (the max- 
imum of the dx2-y2 functional form), but instead at "hot- 
spots" at the intersection of the FS with the antiferro- 
magnetic Brillouin zone (BZ) boundary [18]. This, along 
with the fact that the antiferromagnetic and supercon- 
ducting phases in the highest-Tc samples (x = 0.15) are 
in very close proximity to each other and may in fact 
coexist, [20] led to the conclusion that this high-energy 
pseudogap was due to the effects of antiferromagnetic 
or other orders with characteristic wave vector (tt, tt) 
[18]. As the signatures of boson-electron interactions 
are somewhat generic, it may be difficult to distinguish 
interactions with different types of modes (for instance 
phonons vs spin fluctuations) if the regions of the BZ, 
where such effects are observed to be the strongest, over- 
lap. In the context of Eq. 1, this raises the riuestion 
of the fc-independence of the self-energy implied in this 
type of analysis, especially in the n-type materials where 
the momentum-dependent scattering appears to be very 
strong even near optimal doping [18]. 

Keeping the above considerations in mind, we present 
an ARPES line-shape analysis of the electron doped (n- 
type) cuprate superconductor Ndi.85Ceo.i5Cu04. We use 
essentially same methodologies that have been employed 
in the study of the j>type compounds as we believe that 
a direct comparison has utility despite the fact that there 
may be possible complications. We find that there is no 
kink along the zone diagonal, in contrast to the p-type 
materials where this is the region of the BZ where it can 
be most unambiguously identified. In the (tt, 0) to (tt, tt) 
direction, a very weak kink in the quasiparticle disper- 
sion can be found. However, it is not straightforward to 



identify this effect with phonons. since competing effects 
cannot be separated unambiguously in this region of the 
BZ. Despite an apparent lack of a mass renormalization 
along the zone diagonal, a very slight drop in the scat- 
tering rate is observed in this direction. No such drop in 
the scattering rate is observed near the (tt, 0) to (tt, tt) FS 
crossing. As it is reasonable that this class of many-body 
effects along the nodal direction in the p-type compounds 
can be identified with effects of phonons, the lack of a 
kink and weak drop in the scattering rate along the same 
direction indicates a weaker electron-phonon interaction 
in Ndi.85Ceo.i5Cu04. 

II. EXPERIMENT 

Angle-resolved photoemission data were recorded at 
Beam line 10.0.1.1 of the Advanced Light Source with a 

photon energy of 55 eV in a glancing-incidence geometry. 
The polarization angles were 45° to the Cu-0 bonds for 
the F to (tTjTt) cut and along the Cu-0 bonds for the 
(tt, 0) to (tt.tt) cut. The energy resolution was typically 
18 meV and angular resolution was 0.14° (corresponding 
to w 0.5% of the Brillouin zone size). This configuration 
is identical to that used in an extensive study of the kink 
in the hole doped compounds [10]. The samples were 
cleaved in situ at low temperatures in vacuum better 
than 4 x 10""'^^ Torr. As there are only minor changes 
to the spectra when entering the superconducting state 
[19], all displayed spectra were taken at low temperature 
(w 15-25 K). No signs of surface aging were seen for the 
duration of the experiment (~ 24 h) . In addition to the 
displayed data, extensive measurements were performed 
on Beam line 5-4 of the Stanford Synchrotron Radia- 
tion Laboratory; results were obtained that are consis- 
tent with those displayed here. 

Single crystals of Ndi.85Cco.i5Cu04 were grown by the 
traveling-solvent grown method in 4 atm. of O2. The 
as-grown crystals are not superconducting and must be 
annealed to remove excess apical oxygen. The displayed 
data come from samples that were heated in flowing ar- 
gon gas for 20 h at 920° C, and then a further 20 h at 
500°C under pure O2 [21]. Resistivity and magnetic sus- 
ceptibility measurements show an onset of the supercon- 
ducting transition at 24 K. 

III. RESULTS 

In Fig. 1, we show the raw data in the form of 2D 
intensity plots of the spectral function, along with the 
corresponding EDO's for fc-space cuts, in the F to (tt, tt) 
and the (tt, 0) to (tt, tt) directions. In Fig. 2, we show 
the same data represented in a more traditional fashion 
with EDO's selected at larger angular intervals (0.5°). 
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FIG. 1. Data from high symmetry directions with E~f=55 eV (a) Image of the spectral function along the F to (tt, tt) direction. 

(b) EDO's from (a), (c) Image of the spectral function along the (tt, — tt) to (tt, tt) direction, (d) EDO's from (c). The lower piece 
of the parabola in (c) is slightly distorted from its intrinsic shape as these data were at the edge of the detector multichannel 
plate and were likely affected by astigmatism in the electron lens. 



In the r to (tt. tt) direction, a large broad feature dis- 
perses out of the background, sharpens to a peak at Ep, 
and then disappears. As pointed out previously [18], this 
is behavior not unlike that found in optimally doped p- 
type compounds. In the (7r,0) to (7r,7r) direction, an 
almost parabolic-shaped band centered around (7r,0) is 
observed. A similar sharpening of the feature as it dis- 
perses to Ep is seen here as well, however, the \k\ < Ikp] 
EDO's seem to show additional structure in the form of 
two separate peaks. This aspect has been also noted pre- 
viously [18,19,22]. 




Energy (eV) 

FIG. 2. £^,=55 eV EDO's representing the data of Fig. 1 in 
a more traditional fashion, (a) (0.257r, 0.257r)— >(0.657r, 0.657r). 
(b) (7r,0)^(7r,0.467r). The EDO's are offset for clarity. 

The MDC's from the raw false intensity plots in Fig. 
la and Ic were fitted to Lorentzians as detailed above. 
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FIG. 3. MDO's from Fig. 1 at w = and w = 150 meV 
binding energy along the two main symmetry directions, (a) 
F— »(7r,7r). (b) (tt, 0)— »-(7r, tt). The spectra are offset for clarity. 

Representative fits are shown in Fig. 3 at w = and co = 
150 meV slices. At energies higher than those displayed, 
the MDC line-shape begins to difiier more significantly 
from a simple Lorentzian. 

At the displayed relevant lower energies, the 
Lorentzian functional forms fit the MDC's relatively well, 
showing that, although there may be large momentum 
dependence in the self-energy in other parts of the BZ, 
the changes across the displayed regions are minimal at 
low energy. The fact that the fits begin to break down 
at higher energies indicates that momentum dependent 
effects may become relevant. These fits to the peaks po- 
sitions are displayed in Fig. 4 with the corresponding 
energy widths S" shown in Fig. 5. Horizontal and verti- 
cal errorbars in the Figs. 4 and 5, respectively, represent 
the uncertainties in the values from the Lorentzian fits. 
In Fig 4a we show the dispersion along the zone diagonal 
plotted from the momentum space peak position. As seen 
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by comparison with the hncar fit, there is only a gentle 
slope change and no sign of a sharp mass-renormalization 
at low energy. This is in sharp contrast to the p-type ma- 
terials where this is the region of the BZ that shows the 
most unambiguous effect at « 70 meV. 

The behavior along the (tt, 0)-(7r, tt) direction differs 
from that of the zone diagonal. Here, the smooth disper- 
sion deviates at low energy from the parabolic fit in the 
region indicated by the arrow in Fig 4b. The deviation 
in the dispersion is greater than the uncertainty in the 
fit. The small mass renormalization around 55 meV has 
the same energy scale as the dip structure in the EDC in 
Fig. 2. This is close to the same energy scale as the mass 
renormalization onset in the p-type compounds. If one 
is to identify this apparent mass change with coupling to 
a bosonic mode, a simple analysis for the dimensionless 

coupling parameter A' = ^ — 1 gives a value A' ~ 0.3. 
This may be compared with values » 0.5 extracted via 
the same method from the zone diagonal direction of the 
p-type materials near optimal doping [23]. In this sim- 
ple analysis, we used as a value for Vp the velocity from 
the dispersion above the expected kink energy and the 
VF from a low -w fit, as was done in the analysis for the 
p-type materials [10]. 

Closer to the spirit of how a mass renormalization 
is typically defined would be to use a w ^ extrap- 
olation of a function fitted to the dispersion at higher 
energy (oj > 100 meV). Analyzing the data in this 
fashion gives velocities of "iPp^^ {oj — > 0) = 4.3 x 10^ 
m/sec (2.3 eY-a/hir) for the r-(7r,7r) FS crossing and 
iPpff{uj ^ 0) = 3.4 X 10^ m/sec (1.8 eV-a/h-K) for the 
(tt, 0)-(7r, tt) FS crossing. Using this value for v^p for the 
(tt, 0)-(7r, tt) cut gives a smaller A' of « 0.1. However, 
the values for A' near (tt, 0) arrived at in this fashion are 
not readily comparable to those of the hole doped com- 
pounds along the zone diagonal as their dispersion there 
does not appear to recover to something resembling an 
unrenormalized dispersion within the experimental win- 
dow. This is most likely a consequence of the fact that 
for straightforward electron-boson coupling one expects 
that the dispersion recovers to the bare one within about 
5 times the bosonic band width. In the (tt, 0) region of 
the BZ, the local bare electronic band width may be nar- 
rower than this bosonic bound, so one has to regain a 
parabolic dispersion more quickly. 

We note that finite resolution effects should not cause 
an inaccurate parametrization except very close to the 
Fermi energy. Close to Ep where the Fermi function 
differs significantly from unity, finite energy resolution 
(A_E) weights the signal at each individual uj point dis- 
proportionately towards higher energy so that only fits to 
w higher than /S.E are reliable and intrinsic. The sharp 
change to higher velocity at w < 15 meV in the V to 
(tt, tt) direction is obviously explained by effects of this 
kind. As will be shown later, the data from this energy 
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FIG. 4. Dispersions fitted from the spectral function anal- 
ysis, (a) F— >(7r,7r). The sharp upturn at low energies (w < 15 
meV) is believed to be an artifact of the finite energy reso- 
lution, (b) (tt, 0)^(7r, tt). The smooth lines are simple fits 
(linear and parabolic for (a) and (b) respectively). The hor- 
izontal bars represent the uncertainty in the peak position 
from the Lorcnztian fits. The arrow marks the position of the 
low-energy kink. 

region show also a sharp decrease in scattering rate along 
the zone diagonal. We consider this also to be an artifact 
of finite energy resolution, although the reason such an 
effect does not show up in the (tt, 0) to (tt, tt) direction is 
an open question. 

As detailed above, the widths of the MDC Lorentzians 
can be parametrized as T,"{k,u})/v^. Like the case of A', 
in the absence of a definitive independent measure for Vp, 
we use for i/p-^ ^ the v from a w — > extrapolation of a fit 
to the higher-energy dispersion to get a rough estimate 
for S". For instance, in the (tt, 0)-(7r, tt) direction, we 
use the value for the if^ of the fitted parabola. This fit, 
although displayed as extrapolated to w = 0, does not 
use input data for oj < 100 meV. 
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FIG. 5. Energy widths E" as fitted from spectral function 
MDC's and multiplied by velocity, (a) r^(7r,7r) The solid 
and dashed line are power -law fits [ to S"(a;) = 'S"{0) + Au}"] 
with a = 1.55 and 1, respectively. The arrow marks the point 
of low- energy drop in S". (b) (tt, 0)^(7r, tt). The solid line 
is a linear fit. Tiie vertical bar at 0.18 eV marks the energy 
beyond which the MDC width fits are not reliable. Vertical 
bars represent the uncertainty in the fitted Lorenztian width. 

Along the r-(7r, tt) direction (Fig. 5a) S"(ti^) is best fit 
(solid line) in the higher energy region (0.4 eV> ui >0.09 
eV) to a power law form, S"(a;) = 5]"(0) + Au'^, with 
E"(0) = 0.14 eV and a = 1.55. Note, that the best fit to 
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a linear dependence (dashed line) with a S"(0) of 0.11 eV 
does not fit nearly as well. We observe a low-energy drop 
in E" at approximately 0.07 eV. Although the drop is 
subtle, we believe the effect to be intrinsic as it was seen 
in a number of samples and is outside the uncertainty in 
the fits as represented by the error bars. Note that this is 
close to the same energy that the mass renormalization 
is found at in the (tt, 0)-(7r, tt) direction. The precipitous 
drop in S" for w < 15 meV is found from fitted MDC's 
that gave the artificial upturn in the dispersion in Fig. 
4a and therefore appears to be an artifact of finite energy 
resolution. 

In contrast, along the (tt, G)-(7r, tt) direction, 5]"(cj) is 
best fit by a linear dependence on w with S"(0) = 0.08 
eV with no signature of a low-w drop in the scattering 
rate. Note that due to the much smaller local electronic 
band width in this region of the BZ, the MDC Lorentzian 
fit breaks down at energies higher than « 0.2 eV as may 
have been expected from a visual inspection of Fig. Ic. 
It is interesting to note that S"(0), extracted in this way, 
is very different for the two directions. Also note that the 
finite energy (AE) and momentum resolution (Afc) only 
give a minor broadening contribution to S" as AS" = 

^J{dE^/dk ■ AA?)2 + (A£;)2 « 28 meV. 

IV. DISCUSSION 

As discussed above and in Ref. [10], due to its ubiquity, 
temperature dependence, and energy scale, the strong ve- 
locity renormalizations in the p-type materials has been 
attributed to an interaction of the charge-carriers with 
phononic degrees of freedom. In the n-typc compounds 
we see some evidence of similar effects although impor- 
tant differences exist. As detailed above, there is no kink 
along the zone diagonal, which is the region of the BZ 
that shows the most unambiguous effect in the hole doped 
materials. There does appear to be a kink along the (tt, 0) 
to (tt, tt) direction. However, as we will discuss below, the 
interpretation of this feature is not straightforward. 

If there exists a symmetry reducing order (or associ- 
ated fluctuations) with characteristic wave vector (7r,7r), 
then one may expect to see two features in the EDO's in 
portions of the BZ where the unreconstructed band struc- 
ture is below Ef in momentum regions straddling both 
sides of the antiferromagnetic BZ boundary since the ef- 
fect of such an order is to reflect features across the an- 
tiferromagnetic BZ boundary, i.e., spectral features from 
{h, k) would be expected to show up at {—h + n, ~k + n). 
For instance, EDC's from near the Fermi surface near 
(tt, O.Stt), would have features reflective of the unrecon- 
structed bands well below Ep from near (O.Ttt, 0). In 
contrast, EDC's from the band near (0.457r, 0.457r) will 
not have features that derive from near (O.SStt, O.SStt) as 
this region of the BZ in the unreconstructed band struc- 



tiire is well above Ep. In this way, such a BZ reducing 
order may be expected to greatly modify the spectral 
function near (tt, 0), but not as much near (7r/2,7r/2) for 
an unreconstructed hole-like Fermi surface centered at 
(tt, tt). Two such features in the EDC can give a kink in 
the effective dispersion if an MDC analysis is performed. 
As it has been previously noted that there are two fea- 
tures in the EDC's near the (tt, 0) position [18,19,22] (see 
Fig. 2b) and there is evidence for fluctuations of some 
kind of (tt, tt) symmetry reducing order, we do not be- 
lieve that there is sufficient evidence to unambiguously 
assign the {tt, 0)-(7r, tt) kink feature, found by an MDC 
analysis, the same origin as the kink observed in the p- 
type compounds. We believe there is no means to unam- 
biguously separate the (tt, tt) symmetry reducing effect. 
In this regard, the kink found along (tt, 0)-(7r, tt) in the 
present work appears to be of a different sort than the 
zone diagonal kink found in the p-type. The current case 
may be an example of a mass renormalization due to the 
effects of generic residual antiferromagnetism (i.e. not a 
resonance mode) or other symmetry reducing order with 
characteristic wavevector {tt,tt), and therefore its obser- 
vation does not directly bear on the debate regarding 
the relative merits of the phonon vs magnetic resonance 
mode scenarios for the spectral function of the cuprates. 

Despite the uncertainties in the MDC analysis one can 
make a strong statement about the lack of a kink along 
the zone diagonal. By the same condition as given above, 
in this region of the BZ, we do not expect to see (tt, tt) 
features (nor are any observed) as the unreconstructed 
Fermi sea is only filled on one side of antiferromagnetic 
BZ boundary. Visual inspection of the MDC fitted dis- 
persion in the 0.2 eV > w > eV range shows there is 
no feature resembling a sharp kink in the dispersion [24] . 

The dependence of the scattering rate on frequency 
along the zone diagonal is different from the linear depen- 
dence report(xl in the optimally doped p-type materials 
[25] (note that some groups have found possible devia- 
tions from linear [26]). This stronger u)^-^^ dependence 
in the electron doped compounds compares more closely 
to the quadratic temperature dependence of their resis- 
tivity that to the linear temperature dependence of the 
hole doped materials at optimal doping [27,28]. 

Our results for S" may also be compared to the lin- 
ear dependence of the low-energy (>1000 cm~^) scat- 
tering rate that has been found by an inversion of the 
optical conductivity [29,30]. The conventional wisdom is 
that zone diagonal states give the primary contribution 
to transport due to their supposedly higher vp. At first 
glance this is not consistent with our finding here. How- 
ever, as we have found for the electron doped material 
that the vp at the r-(7r,7r) FS crossing is not greatly 
different from that found at the (tt, 0)-(7r, tt) FS crossing 
[this is another consequence of the fact that the (tt, 0) 
saddle point is « 300 meV below Ep], the (tt, 0) re- 
gion (where a linear dependence is found) may contribute 
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somewhat equally to the optical scattering rate. More- 
over, as may be seen in Fig. 5a, the scattering rate along 
the zone diagonal can be fit approximately by a linear 
dependence above 130 meV and it may be that it is this 
energy which is contributing the most to the optical ex- 
citations for E-y >1000 cm~^ (w 100 meV). Optical ex- 
periments also find a low-energy drop in the scattering 
rate below an energy (650 cm « 70 meV) very close 
to that where we find the depression in S" and a mass 
renormalization [29,30]. 

The energy at which we observe the mass renormaliza- 
tion along the (tt, 0) to (tt, tt) direction is the same energy, 
to within experimental uncertainty, that we observe the 
drop in the scattering rate along F - (7r,7r). Although 
these energy scales are consistent with each other, this 
observation is inconsistent with Kramers-Kronig consid- 
erations. Since S' and S" are real and imaginary com- 
ponents derived from the same response function, they 
should be related via a Kramers-Kronig transform due 
to causality. If T,' along a particular cut shows a kink, 
then S" should show a corresponding drop. The rea- 
son why such a behavior is not observed is unclear, but 
may be related to momentum dependence in the self- 
energy masking such effects. If the self-energy is mo- 
mentum dependent, the dispersion and width obtained 
from an MDC analysis cannot be simply related to S' 
and S", respectively. It does seem that the qualitative 
conclusion that there is little renormalization along the 
zone diagonal remains valid. The situation for the (tt, 0)- 
(tt, tt) direction is less clear, as the mass renormalization 
is present here and the situation to compare to in the 
p-types near (tt, 0) more uncertain. 

As detailed above, it seems reasonable to associate the 
mass renormalization on the p-type side of the phase 
diagram with a strong coupling of charge carriers to 
phononic degrees of freedom. Specifically, the principal 
phonon has been conjectured [10] to be the oxygen half- 
breathing mode, which is seen to have anomalous soft- 
ening and spectral weight changes in the relevant energy 
range near the doping induced metal-insulator transition 
as probed by neutron scattering [31 34]. A nmnber of 
recent studies have found similar signatures of phononic 
anomalies in the electron doped materials [35,36]. Kang 
et al. found changes with doping in the generalized 
phonon density of states around « 70 meV by neutron 
scattering. d'Astuto et al. measured phonon dispersions 
via inelastic x-ray scattering and assigned the softening 
they found in the 55 - 75 meV energy range to the same 
oxygen half-breathing mode that anomalies are found in 
in the p-type materials. These studies give evidence for 
phononic effects in the electron doped materials that are 
somewhat similar to those found in the hole doped com- 
pounds. 

There are a number of differences that do exist in 
these anomalies that may provide a starting point for 
understanding why similar signatures in the boson spec- 



trum may not cause similar ecfccts to appear in ARPES. 
Although the biggest changes in the phonon density of 
states probed by Kang et al. happen at similar doping 
levels in La2-2;Sr2;Cu04 and Nd2-2;Ce2;Cu04 (x ?a 0.04), 
the doping levels are at very different relative positions in 
the phase diagram, with x = 0.04 being still well into the 
antifcrromagnctic (and possibly more insulating) phase 
for the electron doped compound. As such modifications 
in the phonon spectrum may be associated generally with 
screening changes (and hence electron-phonon coupling) 
with the onset of metallicity, this demonstrates the pos- 
sibility that the changes in the NCCO phonon spectrum, 
although superficially similar in the electron and hole 
doped materials, are in some sense different. d'Astuto 
et al. found that although the softening of the phonon 
dispersion appears in a somewhat similar way at a simi- 
lar energy scale as in the p-types compounds, differences 
do exist in the shape of the anomaly in the phonon dis- 
persion. Moreover, on general grounds, since the pur- 
ported soft phonon is the oxygen half-breathing mode, 
one may naively expect a weaker coupling for this mode 
with electron doping, as Madelung potential considera- 
tions [37] indicate that doped electrons will preferentially 
sit on the Cu site, whereas doped holes have primarily 
oxygen character. 

Kang et al. have mentioned [35] that it would be inter- 
esting to look for mass renormalizations in the ARPES 
spectra of the n-type underdoped compoimds that they 
first detect phonon anomalies in. However, we believe 
that if the electron-phonon coupling is of the same na- 
ture in the n— and p— type materials it should reveal 
itself in the ARPES spectra of the highest-Tc electron 
doped samples, as it does in the hole doped materials. 
Moreover, the issue of looking for mass renormalizations 
and at peak widths may not even be relevant at the ex- 
tremely low dopings (x = 0.04) where Kang et al. see the 
largest changes in the phonon density of states, as there 
are only regions of the BZ with enhanced near-_Ei? spec- 
tral weight and no well-defined electronic peaks [38]. The 
qiicstion of what the electron-phonon coupling parame- 
ter is as one enters the Mott insulating state at very low 
doping, may not even be particularly well posed, as the 
near- i^i? spectral weight vanishes as x = is approached. 

V. CONCLUSION 

A modern spectral function analysis reveals very dif- 
ferent and likely weaker signatures of possible bosonic 
effects in the ARPES spectra of Ndi.85Ceo.i5Cu04 along 
the two main symmetry directions of the BZ as compared 
to the hole doped compounds. However, unlike the p- 
type compounds, we believe that one cannot distinguish 
the various contributions to the self-energy into separate 
effects, as the simple analysis yields results which are not 
consistent with Kramers-Kronig considerations. Despite 
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the complexities inherent in assigning various features in 
the spectra to distinct sources, a direct comparison be- 
tween the p- and n-type compounds does have utiUty, as 
we have shown the lack of a mass rcnormalization along 
the zone diagonal. Our analysis shows that whatever 
is causing the distinct effect in the hole doped materi- 
als manifests itself less strongly on the n-type side, per- 
haps due to a weaker electron-phonon coupling or mask- 
ing effects from other interactions. Although such ob- 
servations arc consistent with a weaker electron-phonon 
coupling in the electron doped cuprates, a strong dis- 
crepancy exists between these measurements and a num- 
ber of scattering studies that point to strong electron- 
phonon coupling. Perhaps, some of the intrinsic differ- 
ences between electron and hole doped materials may 
provide a route towards reconciling these different mea- 
surements and thereby giving insight into the underlying 
phenomenon. 
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